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In single particle coherent x-ray diffraction imaging experiments, performed at 
x-ray free-electron lasers (XFELs), samples are exposed to intense x-ray pulses to 
obtain single-shot diffraction patterns. The high intensity induces electronic dy- 
namics on the femtosecond time scale in the system, which can reduce the contrast 
of the obtained diffraction patterns and adds an isotropic background. We quan- 
tify the degradation of the diffraction pattern from ultrafast electronic damage by 
performing simulations on a biological sample exposed to x-ray pulses with different 
parameters. We find that the contrast is substantially reduced and the background is 
considerably strong only if almost all electrons are removed from their parent atoms. 
This happens at fluences of at least one order of magnitude larger than provided at 
currently available XFEL sources. 

Modern crystallography has developed efficient tools for protein structure determination 
up to atomic resolution provided that sufficiently large high-quality crystals are available 
[TJ|2]. Presently there are two major limitations for the future development of this approach, 
namely crystallization and radiation damage. It is well known that the majority of protein 
macromolecules, especially membrane proteins, do not crystalize. At the same time, radia- 
tion damage at conventional x-ray sources limits the available resolution of non-crystalline 
biological samples to few tens of nanometers [3]. 

To break this radiation barrier limit in the structure determination of non-crystalline 
samples it was suggested [I] to use high power x-ray free-electron lasers (XFELs) [5H?] in 
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so-called 'single-molecule' experiments [HI [9]. In this technique, single biological particles are 
injected in the strong coherent XFEL beam in a random orientation, and their diffraction 
patterns in the far field are measured. By collecting a large series of such diffraction patterns 
from a reproducible sample and applying orientation determination techniques [101412"] . the 
full three-dimensional (3D) diffraction pattern in reciprocal space can be assembled. Using 
phase-retrieval algorithms [TBI [H] a 3D image of the electron density of the particle can be 
reconstructed in the final step. With XFEL pulses shorter than few tens of femtoseconds, 
no radiation damage of the sample due to Coulomb explosion is expected during the pulse 
propagation [4J. In this case the diffraction pattern should contain, at least in principle, 
information about the sample in its initial undamaged state. This method has been success- 
fully applied recently to protein nanocrystals [15] and viruses [16] with sizes of few hundred 
nanometers. However, the resolution of non- crystalline samples was found to be far below 
the expected sub-nanometer limit. 

To increase the resolution and to allow the imaging of smaller samples, ideally individual 
proteins, the photon fluence should be increased by several orders of magnitude. In this high- 
power regime, however, atoms are known to undergo strong ionization [T7J HE]; we can thus 
expect that the electronic density will be altered considerably during the scattering process. 
This raises the question how much the diffraction pattern is affected by this electronic 
damage, and if there is an upper limit to the fluence, beyond which the quality of the 
diffraction pattern decreases substantially and prevents a successful reconstruction. This 
problem was first addressed in Ref. [19], which was focused on adapting phase-retrieval 
algorithms to incorporate the effect of electronic damage. In this letter, we quantitatively 
analyse the degradation of the diffraction patterns for a wide range of incoming x-ray pulse 
parameters. 

In a 'single-particle' experiment, the time-integrated intensity at the detector can be 
expressed as 



where the scattering amplitude per unit flux is given in the kinematical approximation by 
A(q, t) = fi(q, t) exp(«qRi(£)). Here, q is the scattering vector, j(t) is the intensity of the 
incoming pulse, and fi(q,t), R«(t) are the time-dependent form factor and position vector 
of the i-th atom, respectively. In the following we approximate the atomic form factors to 
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be spherically symmetric and consider only small samples with a size of tens of nanometers. 
We also assume a homogenous and coherent illumination over the sample area, and neglect 
retardation effects [20] in Eq. ([I]). Note that the recorded diffraction pattern described by 
Eq. has the form of an incoherent sum over instantaneous coherent diffraction patterns 
|A(q, t)\ 2 . Clearly, under certain experimental conditions, this summation can decrease the 
contrast of the measured diffraction pattern. 

The radiation damage induced by intense XFEL pulses is basically a two-step process. 
In a first step, the sample is strongly ionized on a time scale of few femtoseconds [JRJ Elh 
[21]. The ejected electrons leave the sample, which thereby accumulates a space charge 
[25] . In a second step, the system disintegrates through a Coulomb explosion within tens 
to hundreds of femtoseconds (U EJJ- This results in a motion of the atoms from their 
equilibrium positions, Rj = Rj(t), which seriously degrades the diffraction pattern through 
the incoherent summation in Eq. ([T| . The use of pulses with a duration of few femtoseconds 
should prevent the degradation from the Coulomb explosion [3], but the electronic damage 
cannot be avoided. Obviously, the ionization lowers the atomic form factors, which reduces 
the scattered signal and thus the achievable resolution. In addition, it introduces an explicit 
time-dependence in the form factors, which can also reduce the contrast of the diffraction 
pattern. 

For our analysis, we assume that a large number of single-shot patterns have been mea- 
sured, aligned, and averaged perfectly. Since we are interested only in electronic damage, 
we also require that the atomic positions do not change during the pulse propagation, and 
that the sample is perfectly reproducible for all shots. With these assumptions, an average 
of Eq. (jl]) over many pulses gives 



Here, Iq = J j{t)dt is the total fluence of the x-ray pulse, the brackets (. . .) denote averaging 
over many pulses and the summation is performed over all atoms in the sample. A normalized 
form factor matrix F with matrix elements 



is introduced in Eq. ([2]). We assume here that all FEL pulses are identical |26j, so that the 

intensity j(t) can be taken out of the average in Eq. pi). Note that the time integration 
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in Eq. Q that causes the degradation of contrast is now encapsulated in the form factor 
matrix F. 

To proceed with the analysis, the time- dependent form factors can be written as a sum 
fi(q,t) = (fi(q,t)) +5fi(q,t) with average values (fi(q,t)) and pulse-to-pulse fluctuations 
Sfi(q,t) (by definition (Sfi(q,t)) = 0). We can then express the average of the form fac- 
tor product in Eq. @ as (#(?,*)/;(?,*)) = (f*(q, t)) (f^q, t)) + (Sft(q,t)Sfj{q,t)). If we 
finally assume that the ionization of atoms at different positions i, j is uncorrelated [21 J, the 
second term reduces to (5f*(q,t)Sfj(q,t)) = Sij{\5fi(q,t)\ 2 ). The matrix F then naturally 
decomposes into two contributions, 

F i3 {q) = W ij (q) + 6 lJ B i (q) , (4) 
where W(g) is an Hermitian matrix [27] 

WiM) = jj iW(/*(g,i))(/i(9,i))dt , (5) 
and the vector B(g) is defined as 

B l (q) = jjj(t)(\6f t (q,t)\ 2 )dt. (6) 

Substituting Eq. Q in Eq. ^ we obtain for the average intensity 

(/(q)) = / w (q) + / B (g) , (7) 

where Iw(q) = h Wij ex P[ ?c l(Rj ~~ an d ^B(q) = Iq J2i B%{q). Our derivation shows 
that the average intensity (/(q)) in a 'single-particle' experiment can be presented as a sum 
of two contributions. The structural term, J\y(q), is determined by the average form factor 
values (fi(q,t)), and the background term, /b(<?), by their fluctuations (\5fi(q,t)\ 2 ). 

The structural term Iw(q) i n Eq. (0) can be presented as a sum of coherent modes [T9"] . 
To show that, we formally decompose the matrix W(g) as [28] 

W tJ (q) = J2 c ^ v * a> (M a) ( ( l) > (8) 

a 

where c a , v^ a ^ are the eigenvalues and normalized eigenvectors of W(g). With this decom- 
position, we have for I\y(q) in Eq. ^ 

/w(q)=/fl^c tt ( ? )|A( a )(q)| 2 , (9) 




FIG. 1: (Color online) Averaged two-dimensional diffraction patterns of the biological sample in 
logscale for an x-ray pulse of 5 fs FWHM and 3.1 keV photon energy, (a) Diffraction pattern of 
the neutral sample without electronic damage at a fluence of 10 14 photons//im 2 . (b, c) Diffraction 
patterns of the ionized sample at a fluence of 10 photons//im 2 (b) and 10 16 photons/^um 2 (c). 
(d) Intensity profiles for (a-c) along the black lines. The intensities in (d) have been rescaled to be 
equal at q = 0. 

where A^(q) = £\ v i f\q) exp(zqR, ; ). An inspection of Eq. ^ shows that the structural 
contribution i\y(q) to the averaged intensity (/(q)) is described by a sum of coherent modes 
I -A (q)| 2 - Each mode produces a diffraction pattern with full contrast with the weight 
given by the eigenvalue c a (q). This summation reduces the contrast of the final diffraction 
pattern depending on the incoming pulse parameters. The magnitude of the contribution of 
the modes can be described by the parameter [29] 

_ Tr(W(g)W(g)t) Z a <? a (q) 
m (TrW(g)) 2 ~(£ Q c Q (g)) 2 ' {W) 
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where TrW = J2i Wu denotes the trace operation. By definition, < ((q) < 1, and ((q) is 
equal to one only if there is a single mode at a certain value of q. 

The background term /b(<?) in Eq. (0 gives rise to an additional isotropic background. As 
a result of our assumption that the ionization of atoms at different positions is uncorrelated, 
this background does not contain any structural information and cannot be included in the 
mode decomposition described by Eq. To quantify the relative contribution of this 

background we introduce the ratio 

between the background and the angular average of the structural term Iw(q) i n Eq- 
Here, the brackets (. . .)^ denote angular averaging of the diffraction pattern. 

To analyze the effect of the electronic damage, we have simulated [30J averaged diffraction 
patterns using Eq. ([7]) for a human adenovirus penton base chimera [31]. The sample has 
a dodecahedral shape with a diameter of 27 nm and contains about 200,000 non-hydrogen 
atoms, giving a mass density of about 0.5 g/cm 3 . In our simulations we considered an incom- 
ing x-ray beam with a photon energy of 3.1 keV (corresponding to the wavelength A = 4 A), 
and 12.4 keV (A = 1 A) focused to a 100x100 nm 2 area at the sample position. The detector 
was positioned 100 mm downstream from the sample and had a size of 200x200 mm 2 for the 
soft, and 400x400 mm 2 for the hard x-ray energies. This provided a maximum scattering 
angle of 45° and available resolution d max =5.2 A for the soft, and 63.4° (d max =0.95 A) for 
the hard x-ray energies. 

The averages of the form factors (fi(q,t)) and their fluctuations (\Sfi(q,t)\ 2 ) in Eqs. ^ 
and Q were evaluated with a rate equation approach [22]. We included ionization from 
photoionization, Auger decay and electron impact ionization into the rate equations [32] ■ In 
our simulations, only carbon, nitrogen and oxygen atoms were considered. The contribution 
from sulphur and hydrogen atoms can be safely neglected; the former appear only in trace 
quantities, while the latter hardly scatter. 

Simulated diffraction patterns for an x-ray pulse with 5 fs FWHM, 3.1 keV photon energy 
and fluences of 10 14 and 10 16 photons//xm 2 are shown in Fig. [TJ We can immediately see 
that the electronic damage reduces the contrast of the diffraction pattern in comparison 
with the undamaged sample (Fig. [l](a)). For the lower fluence (Fig. [l]^b)), this mainly leads 
to a slight smearing of the deep minima of the diffraction pattern. At the same time, for 
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FIG. 2: (Color online) (a, c) Contrast C((7max)> and (b, d) relative background contribution r(g max ) 
calculated at (/max corresponding to a resolution of 5.2 A at 3.1 keV photon energy (a, b) and 
a resolution of 0.82 A at 12.4 keV photon energy (c, d). The labels indicate the approximate 
fluences for the saturation of the ionization of the respective shells: Is of neutral carbon, 2s of 
C 2+ (ls°2s 2 2p 2 ), and 2p of C 4+ (ls°2s°2p 2 ). The dashed line in all figures corresponds to the 
fluence of 10 14 photons/^m 2 presently achievable at XFELs. 

the higher fluence (Fig. |1Jc)), most of the details of the diffraction pattern are smeared out. 
To determine the relative contribution of the multimode decomposition and the back- 



ground to the diffraction pattern, we have calculated the parameters ((q) and T(q) [Eqs. (10) 



and (11)] at the momentum transfer value g max = 2n/d max for a wide range of x-ray pulse 
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parameters (see Fig. [2]). At the soft x-ray energy (Fig. |2[a)), C(<7max) is close to one up to 
the fluences of 10 15 photons//im 2 when all electrons are photoionized. At higher ffuences 
it rapidly drops down to a value of 0.9. At hard x-ray energies (Fig. [2^c)), the contrast 
drops in two steps. First, slightly (down to a value of 0.95), at the fluences from 10 15 to 
10 16 photons/yum 2 when core electrons are efficiently ionized. Then, strongly (to a value of 
0.8), at extreme fluences higher than 10 17 photons/ /im 2 when all electrons are removed from 
the atoms. It follows from our simulations that the contrast is substantially reduced only 
if almost all scattering electrons are removed from their parent atoms. This will not cause 
any problems with the presently available fluences but will become a substantial factor at 
fluences at least one order of magnitude higher. 

The dependence of the background contribution r(g max ) on the pulse parameters 
(Figs. [2^b,d)) is qualitatively similar for the soft and hard x-ray energies. The background is 
negligible for the lowest fluences and rises continuously with increasing fluence. It becomes 
significant (more than 10%) at fluences higher than 10 14 (10 15 ) photons//zm 2 at the soft 
(hard) x-ray energies. At high fluences, r(g max ) ranges from 0.4 to 0.9 at soft x-ray energies 
and from 0.5 to 2 for hard x-ray energies for pulse durations from 1 fs to 40 fs. For hard 
x-ray energies and fluences of about 10 16 photons//im 2 , the background is reduced for pulses 
shorter than 3 fs. We attribute this to formation of hollow atoms |22j, which reduces the 
form factor fluctuations due to the K-shell depletion. 

As an important outcome of our simulations we could estimate an effect of electronic 
damage on the scattered signal with increased fluence (see Fig. [3|a,c)). For a soft x-ray 
pulse with 10 14 photons//im 2 the sample scatters one order of magnitude less photons than 
a corresponding undamaged sample. In contrast, for high x-ray energies, both signals are on 
the same level (see Fig. [3^c)). In the first case, the fluence is sufficient to remove most of the 
(core and valence) electrons from the atoms (see Fig. |2^a)), while for high x-ray energies the 
photoionization cross sections are substantially smaller (see Fig. [2^c)). Furthermore, once 
the fluence increases beyond the level of 10 13 (10 15 ) photons//xm 2 for soft (hard) x-rays, 
an increase of the scattered signal by one order of magnitude requires an increase of the 
incoming intensity by three orders of magnitude due to excessive ionization (compare with 
[22]). 

To estimate limitations to the achievable resolution introduced by the electronic damage 
we analyzed the q— dependence of the structural term (Iwil))^ an d the background ratio 
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FIG. 3: (Color online) (a, c) Angular averaged scattered signal {Iw(q)},, and (b, d) relative 
background contribution T(q) for a 5 fs FWHM x-ray pulse at different fluences; (a, b) 3.1 keV, 
(c, d) 12.4 keV photon energy. Dashed curves correspond to the intensity scattered by a sample 
without electronic damage. Dotted lines in (a,c) denote the cut-off of 10~ 2 photons per Shannon 
angle, and the labels denote the achievable resolution. The curves have been smoothed with a 
Gaussian filter to remove high-frequency oscillations from the speckle pattern. 

r(g) (see Fig. [3]). We assumed that the resolution is determined by the criterion of measuring 
10 -2 photons per Shannon angle [TOT ITT] . At hard x-ray energies, a high resolution below 
1 A can be, in principle, achieved with a fluence of 10 15 photons//im 2 (see Fig. [3](c)). For 
these parameters, the background contribution T(q) is about 5%. For higher fluences, how- 
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ever, the background quickly increases; at a fluence of 10 16 photons//im 2 , the background 
signal is about 20% of the scattered signal already at the resolution of 3 A (see Fig. |3](d)). 
Similar behavior can be observed at soft x-ray energies at lower resolution and fluences 
(see Figs. [3^a,b)).This clearly demonstrates that electronic damage during propagation of 
ultrafast pulses limits the achievable resolution at high fluences. 

In summary, we have shown that electronic damage can significantly reduce the contrast 
of single-particle diffraction patterns and gives rise to an isotropic background if high power 
FEL pulses are used. We have derived quantities to describe these effects, and calculated 
them for a typical biological sample illuminated by an intense x-ray pulse. Our simula- 
tion have shown that at soft x-ray energies sub-nanometer resolution can be, in principle, 
achieved for fluences higher than 10 13 photons//im 2 , however, at fluences higher than 10 14 
photons//zm 2 the background will give a substantial contribution to diffraction patterns. 
For hard x-ray energies, 2.7 A resolution is feasible for presently available XFEL sources. 
No contrast degradation or background contribution is expected in this case. Raising the 
fluence up to the limit of 10 15 photons//zm 2 will allow to reach 1 A resolution limit. How- 
ever, further increase in the fluence will rapidly increase the background contribution to 
diffraction patterns. 

We acknowledge helpful discussions with A. Singer, O. Yefanov, and J. Gulden. Part of 
this work was supported by BMBF Proposal 05K10CHG. 
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Supplement: Implementation of the rate equations 

Here, we present the technical details of the implementation of the rate equations. 

The rate equations are based on the approach developed in [22] • We first define a set of 
states that each atom type can occupy. As such states, we consider the energetically lowest 
states for all different occupations of the electronic shells. For example, carbon can have 
between zero and two electrons in each of the Is, 2s, 2p shells, yielding a total of 27 distinct 
states. We can then calculate the time-dependent occupation probability of the £-th state, 
Pz(t), by solving a set of coupled differential equations of the form 

i>&) = E R tv(t)Pv(t) - Rndt)Pt(t) • (12) 

Here, R^if) denotes the total time-dependent rate of transition from state rj to state £. As 
initial condition, we start with all atoms in the neutral ground state. 
In our simulations, the total rate was composed of four parts, 

R(t) = R photo (t) + R Augcr + R escapc (t) + R trap (t) , (13) 

where RP hoto is the rate of photoionization, R Au s er is the Auger rate, and R esca P e ; Rtrap are 
the rates of secondary ionisation from escaping and trapped electrons, respectively. 

The photoionization rate is calculated from R photo (t) = & photo j(t). The photoionization 
cross sections (5- photo as well as the Auger rates R Au s er are calculated within the Hartree- 
Fock-Slater (HFS) approximation |36j; the explicit equations can be found for example in 
[22]. 

It is known that secondary ionization can significantly change the ionization behaviour 
of the atoms [21], so we also accounted for these effects. We use a different treatment for 
photoelectrons on one hand, and Auger and secondary electrons on the other hand. 

Photoelectrons are assumed to originate in the center of the particle, and leave it instan- 
taneously. We consider the particle as a homogenous sphere. With these approximations, 
the rate of secondary electrons produced by the photoelectrons is given by 

R escapc (i) = 6- impact (^ )- dA ^ phot ° , (14) 

(jit* 

where R, V are the radius and volume of the molecule, and diVphoto/dt is the flux of pro- 
duced photoelectrons. The impact cross sections <7 impac were calculated using the Binary- 
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Encounter-Bethe model [37] with electron orbital parameters obtained from the HFS calcu- 
lations. When a secondary electron is produced, its kinetic energy is set to a constant value 
Eq, which we chose to be 25 eV [2T] . 

Auger electrons and all secondary electrons are assumed to be trapped by the charged 
particle, and to thermalize instantaneously into a homogenous electron gas with a Maxwell- 
Boltzmann distribution of velocities. The temperature of the gas is given by the average 
kinetic energy E of the trapped electrons as E = 3/2kT. The rate of collisional ionization 
is then calculated as 

Here, n trap is the density of trapped electrons, v is the velocity of the particles, and ()y is 
the average with the appropriate Boltzmann factors. 

The electron density and total kinetic energy of the trapped electron gas are calculated 
from the set of equations 

d/7-trap _ \ "> \ r / pAuger I p esca P e /j\ | p tra P 



dt 



(C + (*) + RfZ(t))pUt)Q« (is) 



The index a denotes the atom type (carbon, nitrogen etc.) with corresponding atom density 
g a . The energies E^^ r , E^° l } a and E are the energy of the released Auger electrons, the 
binding energy of the electron that is released through collisional ionization, and the energy 
assigned to secondary electrons produced by photoelectrons, respectively. We should note 
that this model is a simplified version of the model used by Hau-Riege et al. [21], and 
was chosen for its simplicity. However, it is able to reproduce most features of the more 
sophisticated simulations in [21], giving some confidence that it does capture the essential 
physics. 



Solving the differential equations (12) yields the time-dependent occupation probabilities 
P£i(t) for each state & of atom i. Furthermore, from the HFS calculations, we can also obtain 
the form factors /^(g) for each state. Using these, we can calculate the average form factors 
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and their fluctuations as 

</i(?,f)>=£p&(*)/&(?) (17) 
(\SMq,t)\ 2 ) = (\Mq,t)\ 2 ) - U(q,t))\ 2 = I /&(?)! 2 ~ |</i(?,*)>| 2 ■ (18) 



